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We analyze the finite temperature structure of the standard model coupled to gravity with one and two dimen- 
sions compactified (on a circle and a torus). We find that finite temperature effects wash out any vacua which 
exist at zero temperature. We also discuss the possibility of transitions between the four-dimensional universe 
and the lower-dimensional spacetimes. 
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I. INTRODUCTION 

It has recently been demonstrated that the standard model 
coupled to gravity, besides the unique four-dimensional vac- 
uum, may also have vacua stabilized by Casimir energies of 
standard model particles if one spatial dimension is compact- 
ified on a circle |l| or if two spatial dimensions are compact- 
ified on a torus [2]]. The vacuum of the low-energy effective 
theory for the one-dimensional compactification can be either 
anti de Sitter, Minkowski, or de Sitter, whereas for the two- 
dimensional compactification it is necessarily anti de Sitter. 
Such vacua exist for a wide range of neutrino masses. 

In this paper we extend the results of ^ and 101 and explore 
the stabilization of compact dimensions at finite temperature. 
We calculate the free energy and use it to derive formulas 
for the energy density and pressure. We then use Einstein's 
equations to analyze the stability of the compact space. We 
find that increasing temperature washes out any existing sta- 
ble points. 

At a given temperature, only particles of masses smaller 
than this temperature have a considerable finite temperature 
effect on Einstein's equations. This fact, along with the known 
standard model spectrum, allows us to analyze the stability of 
compact dimensions up to temperatures around 100 GeV ~ 
IQis K. 

Finally, we explore the possibility of our universe spon- 
taneously compactifying to a lower-dimensional space time. 
The tunneling rate for producing a bubble of such a com- 
pact geometry can be estimated by calculating the Coleman- 
de Luccia-like instanton |3]. Not surprisingly, we find that the 
rate for such a process is exceedingly small. 



n. COMPACTIFIED SPACETIME AT FINITE 
TEMPERATURE 

In this section, we discuss the stability of compactified 
spacetimes for nonzero temperatures. The action of the stan- 
dard model coupled to gravity is given by. 



action with respect to the metric g^^ yields Einstein's equa- 
tions. 
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where 5(4), Mpi ~ 2.4 x 10^** GcV, TZ, and Crcm are the 
determinant of the 4D metric, Planck mass, Ricci scalar, and 
standard model Lagrangian, respectively. We implicitly in- 
clude the 4D cosmological constant in £fioid- Variation of this 
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where G^i, is the Einstein tensor and is the stress-energy 
tensor containing the cosmological constant term. Starting 
from these Einstein's equations, we investigate the structure 
of the standard model for one- and two-dimensional compact- 
ifications. 



A. One-dimensional compactification 

Here we discuss the finite temperature effects for the case 
of one dimension compactified on a circle. Working in four 
dimensions, we consider the spacetime interval, 

ds^ = 5(3)^^ dxf'dx" + R{tf dcj)^ , (3) 

where 3(3)^,^ is the metric on the noncompact 3D spacetime 
with /i, 1/ = 0, 1, 2, and the compact coordinate (/) € [0, 2tt). 
We assume a homogeneous and isotropic 3D metric. 
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where a{t) is the scale factor of the noncompact space, and 
K = 1, 0, or — 1 corresponds to a closed, flat, or open metric, 
respectively. We choose a stress-energy tensor consistent with 
the symmetries of the metric. 
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where is the energy density, and , pji are the pressures 
in the noncompact and compact space, respectively, all being 
functions of R and p = l/{kBT), with T - temperature. Ein- 
stein's equations (|2|i take the form. 
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FIG. 1 : The RHS of equation (9)i ^3 (^). ^ ^ function of the compactification radius R for normal hierarchy Dirac neutrinos with the lightest neutrino mass 
of 6 x 10-12 GoV (red), 7 x IQ-^^ GeV (orange), 8 X lO^i^ GoV (green), and 9 X IQ-^^ GeV (blue) for temperatures: (a) T = 0, (b) T = 30 K, 
(c) T = 40 K, and (d) T = 100 K, respectively. Note that the ranges in figures (c) and (d) are larger than for the first two plots. Alpi has been set to 1 for 
simplicity. 



Here the dot denotes the derivative with respect to time. This 
yields the equation of motion for R, 
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The left hand side of (|9|l can be written as — □ logi?, where 
the d' Alembertian □ = I?,, 13^. The energy density and pres- 
sures are calculated from the free energy using standard ther- 
modynamic relations i4j]. 
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where V2 is the volume of the noncompact space. The total 
free energy is obtained by summing the contributions from all 
standard model particles' and the cosmological constant term. 



It is given by. 
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where A is the cosmological constant, and f!^^^ is the free 
energy for a particle of mass to. It is obtained by following 
the steps outlined in isj]. 
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where the upper sign corresponds to bosons, the lower one 
to fermions, and d is the number of degrees of freedom. The 
quantity V2p'fil^f^^-^ is the zero temperature Casimir energy and 
V2f^^^ is the finite temperature contribution to free energy.^ 



One has to be careful when including contributions from quarks and glu- 
ons. At temperatures of the order of the meson masses, for example, meson 
fields must be treated as fundamental. 



^ Note that while the energy density ps and pressure pa are mass dimension 
four, the Casimir energy density PQ^^s^aj and free energy density f^,^j are 
mass dimension three. 
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It is straightforward to show that formula (ITjt is free of 
divergences. The 3D Casimir energy density for a particle of 
mass m is given by, 
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where Kn{x) is the modified Bessel function of the second 
kind, V2p'(l[^i^^-^^^^ is the divergent part of the Casimir en- 
ergy, and V2p|^(^{(3-) is the finite regularized contribution. In 
the massless limit formula ( fTsT i reduces to. 
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Note here that. 
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where the superscript "obs" indicates the observed value (i.e., 
A°''" ~ 3.1 X 10"*^ GcV^ lH), while "q.corr." denotes the 
quantum correction in flat space given by. 
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Since Pca^(3)(jiv ~ — 27ri? A'^ '^"'"''-, all the divergences cancel 
and the free energy can be expressed in terms of finite quanti- 
ties using equation ( fT3T l, 



2tt RK 



obs 



particles 
PCasA(3) 



/(3) 



(19) 



The first term in the second line of (fT9] l is the Casmir energy 
including the cosmological constant term and has no temper- 
ature dependence. The second term is the finite temperature 
contribution. 

In summary, the dynamics of the field R is governed by 
equation (|9l), where P3,P3,PB. are given by formulas (fTOt - 
(fT2] |. with the free energy expressed through equation (IT9^ . 
It is easy to check that equation ^ at zero temperature be- 
comes the equation of motion resulting from the dimension- 
ally reduced action, derived along the lines of Jll], and takes 
the form. 
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where the prime indicates the derivative with respect to R. 

It has been shown IH that depending on the neutrino 
masses, their type (Dirac or Majorana), and choice of hier- 
archy, there might exist a 3D anti de Sitter, Minkowski, or de 
Sitter vacuum of the low-energy effective theory. Including 
finite temperature effects has similar implications for all those 
cases, therefore, for illustrative purposes, we discuss only the 



case of Dirac neutrinos with normal hierarchy (figure 1). Note 
that the right hand side of equation (|9]l, i.e., W^{R), is anal- 
ogous to the derivative of the potential, so it should be equal 
to zero at the stationary points. In addition, dW^/dR < at 
the stationary point implies stability, whereas dW^ / dR > 
corresponds to an unstable stationary point. 

At temperatures T <C me — 6 x 10^ K only the photon, 
graviton, and neutrino contributions are relevant in equation 
(fT9] l. Figure 1 (a) shows the plot of WsiR) at T = for sev- 
eral lightest neutrino masses. In the cases for which stationary 
points exist, the one appearing at smaller R is stable, whereas 
the one for larger R is unstable. Those points match the ones 
found along the lines of We checked that for low tempera- 
tures (T < 10 K) the finite temperature effects are negligible. 
However, as soon as the temperature reaches approximately 
20 K, the curves slowly depart from the zero temperature re- 
sult and move upwards. Figure 1 (b) shows the plot of W3 (R) 
for T = 30 K. 

With temperature further increasing, each curve moves up 
and eventually fails to cross zero (figure 1 (c), (d)). At 
very high temperatures, the contributions from other standard 
model particles become relevant. We have investigated the be- 
havior of W3 (R) up to a temperature T = 100 GcV ~ 10^'"' K 
and confirmed that no zero points appear. 

The behavior of W^lR) at high temperatures can be un- 
derstood as follows. At a fixed temperature, only particles 
with Kaluza-Klein mass = iri^ + ri^/R^ < contribute 
to /(3), since the particles with larger masses are Boltzmann 
suppressed. In addition, particles with masses smaller than 
T can be treated as radiation. Then, the finite temperature 
contributions to W3 (i?) coming from and (say pT^ and 
p|f , respectively) cancel because pT^ ~ p|^/2 at leading order. 
Equation (|9]l now becomes. 
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The first two terms on the right hand side are the zero temper- 
ature contributions. The third term is temperature dependent 
and can be estimated as /|^) - {TR)'T^ ^ -T"', where the 
factor {TR) comes from summing over Kaluza-Klein modes. 
Therefore, in the high temperature limit, equation (|9) takes the 
form, 

R^ ( /5casA(3) + (J3RT^^ ' 
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where 0-3 encodes the number of standard model particles 
whose mass is smaller than T. For small R we have 
PcasA(3) ~ ^/R^, SO the first term on the right hand side 
of equation ( |22] ) is negligible for T ^ 1/R. Since for all 
experimentally allowed neutrino masses the zero temperature 
vacuum occurs at i?vac ~ 10^° GeV^^, the high temperature 
limit applies for T ^ 1000 K. This confirms that finite tem- 
perature effects remove any zero points of Wj,{R) and, con- 
sequently, make the vacuum disappear. It also explains why 
the universe would have expanded rather than settle in a three- 
dimensional vacuum existing for small R at zero temperature. 
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FIG. 2: The RHS of equation )36t . W2{b), as a function of the volume modulus b for normal hierarchy Dirac neutrinos with the lightest neutrino mass of 
1.0 X 10-" GeV (red), 1.5 X 10"" GeV (orange), 2.0 x 10"" GeV (green), and 2.5 x 10"" GeV (blue) for temperatures: (a) T = 0, (b) T = 100 K, 
(c) T = 130 K, and (d) T = 200 K, respectively. The ranges in figures (c) and (d) are larger than for the first two plots. A/pi has been set to 1 for simplicity. 



B. Two-dimensional compactiflcation 

Next, we investigate the finite temperature structure of the 
standard model for a toroidal compactiflcation. The four- 
dimensional spacetime interval is given by. 
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where and t = ti + iT2 are the volume and shape moduli, 
respectively. From here on we will assume r 
in all formulas for the reasons discussed in |]2|. We also as- 
sume a homogeneous universe, thus b depends only on time. 
For the noncompact spacetime, we parametrize the Hne ele- 
ment as. 



dS(2) = -dt^ + a{tfdr'^. 
The stress-energy tensor has the form. 
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where p2, P2, and pi, are the energy density, pressure in the 
noncompact space, and pressure in the compact space, respec- 
tively. 

We can now proceed along the same lines as in the previous 
section. The 4D free energy for a particle of mass m is. 
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where Vi is the volume of the ID noncompact space and the 
terms in the last line are the flnite part of the Casimir energy 
density, the divergent piece, and the temperature dependent 
contribution, respectively. As shown in the divergent part 
of the Casimir energy is cancelled by the quantum correction 
to the cosmological constant, just as in the previous case. The 
total free energy density in 2D is therefore given by. 
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In the second Une of equation ( |28] ). /Ocas(2)A is the contribu- 
tion at zero temperature including the cosmological constant 
term, while /(2) is the temperature dependent piece. The en- 
ergy density and pressure are calculated from the total free 
energy as before. 
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and yield. 
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As in the previous case, the right hand side of (136) . i.e., W2{b), 
governs the dynamics of the volume modulus of our compact 
space. Stationary points are given by zeros of W2{b). Similar 
arguments as before show that, in the cases for which station- 
ary points exist, the one appearing at smaller b is stable, and 
the one for larger b is unstable. The T = results match the 
numbers obtained in We checked that those zero points 
also disappear with increasing temperature. As an example, 
figure 2 shows the behavior of (b) for Dirac neutrinos with 
normal hierarchy for different temperatures and masses. 

At high temperature the relation P2 — pj (where 
and P2 are the temperature dependent parts of the energy 



density and pressure, respectively) is satisfied at leading or- 
der One can show at next-to-leading order that P2 — P2 ^ 
{bT^T ~ T^, where the factor (&T)2 again comes from 
summing Kaluza-Klein modes. Consequently, equation 
becomes, 
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where 0-2 is mass dimension three and contains information on 
the number of standard model particles whose mass is smaller 
than T. Thus W2{b) is an increasing function of temperature, 
and for high enough T the zero points are washed out. As 
before, we conclude that there shouldn't have been a "two- 
dimensional" epoch in the history of our universe yet. 



III. TRANSITIONS BETWEEN SPACETIMES 

The possibility of transitions between spacetimes of differ- 
ent dimensionality is certainly interesting to explore [i3n7 Hl0ll . 
In our case, as noted in the previous section, there are no vacua 
at high temperature, thus no transition of this type could have 
occurred in the early universe. This is in accordance with 
observation. Obviously, the current visible universe is still 
four-dimensional, which also agrees with the dynamics of the 
compact dimensions presented in the previous section. 

Nevertheless, it seems interesting to investigate the possi- 
bility of a future tunneling of our 4D de Sitter universe to a 
lower-dimensional anti de Sitter spacetime, which we have 
shown to exist for appropriate neutrino types and masses. 
Since the current temperature is T ~ 2.7 K and the universe 
is cooling down, we can safely calculate the transition rates 
using the zero temperature potential (as was shown in the pre- 
vious section, for T < 10 K finite temperature effects are 
negligible). 



A. 4D 2D transitions 

We first note that a transition of our 4D (false) vacuum to a 
2D X anti de Sitter (true) vacuum (0] is only possible if the 
geometry of two spatial dimensions of our 4D universe is that 
of a "very large" torus. The solutions interpolating between 
4D and 2D x regions are black holes with toroidal horizons 
and such transitions correspond to nucleation of those black 
holes in the 4D de Sitter spacetime. In calculating the tran- 
sition rate we follow the usual method outlined in fl''. The 
prescription is simple: we take any finite Euclidean solution 
connecting one side of the potential barrier with the other, cal- 
culate the instanton action, and use this to estimate the transi- 
tion rate. 



' In the transition rate from a 6D de Sitter vacuum to a 4D X 5^ anti de 
Sitter vacuum was calculated. 
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The equations of motion are derived from relations (l33l - 
( |35T l in the Umit T — > 0. The two independent equations re- 
duce to. 



n^:^_ PCasA(2) 
b ~ (27r6)2Af2 



Pl 



C6 



(38) 
(39) 



where C is a constant determined by initial conditions. Ana- 
lyzing numerically the dynamics of b{t) we confirmed that the 
anti de Sitter vacuum found in [0] is indeed stable. 
The transition rate is given by, 

r ^ rQe-(s.„stM-s<iS4) ^ (40) 
where Fq is a constant, the 4D de Sitter action is fstl, 



SdSi 



-24 TT 



yYobs 



-10 



122 



(41) 



and 5'inst [b] is the instanton action. The equations of motion 
in Euclidean space are obtained by flipping the sign of the po- 
tential in equation ( |38] |, as a consequence of the analytic con- 
tinuation t IT. S'inst [b] can be estimated using a finite nu- 
merical (interpolating) solution of the Euclidean equations of 
motion connecting both sides of the potential barrier (equiv- 
alently, connecting points on both sides of the Euclidean po- 
tential stable point). We numerically found finite solutions for 
the initial conditions &(0) = 6o and = for a range of bo 
values. Various finite solutions correspond to nucleating black 
holes of different mass in the 4D spacetime. One of such so- 
lutions, assuming normal hierarchy Dirac neutrinos with the 
lightest neutrino mass of 10^^^ GeV, is shown in figure 3. 
Using this solution, the Coleman-de Luccia-like instanton ac- 
tion can be calculated as fstl. 
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One can also estimate the transition rate by calculating the 
Hawking-Moss-like instanton action which describes b{t) sit- 
ting at the unstable de Sitter vacuum (at b = bds) as, 

- ^HMinst = WTT^'iMpibds)'' ^ 10" < -SdS, ■ (43) 

The Hawking-Moss-Uke instanton is of the same order as the 
Coleman-de Luccia-like instanton obtained via the interpolat- 
ing solution. Such a result was expected since the interpolat- 
ing solution should be equal to the solution at the top of the de 
Sitter hill in the limit bo bds- In both cases the transition 
rate is. 



(44) 



This is an extremely small number, so it is clear that these 
transitions are highly suppressed. 
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FIG. 3: (a) The solution b(r) of the Euclidean equivalent of equation )38t 
assuming nonnal hierarchy Dirac neutrinos with the lightest neutrino mass 
of 10"" GeV and for the initial conditions 6(0) = 6 X lO^*^ GcV"^ and 
b(0) = 0. (b) The corresponding solution a(r) of the Euclidean version of 
equation (39). 



B. 4D 3D transitions 



We now examine the possibility of transitions from the 4D 
de Sitter (false) vacuum to a SDxS*^ anti de Sitter (true) 
vacuum yj]. In the previous subsection we found that the 
instanton action corresponding to the transition to a lower- 
dimensional universe calculated via the Hawking-Moss-like 
instanton gives a result of the same order as the interpolating 
solution. Here we proceed the other way around and first es- 
timate the transition rate to a 3D spacetime by calculating the 
Hawking-Moss-Uke instanton. Equations of motion at zero 
tempereture can be obtained from relations (|6)-(l8) by taking 
T — s- as in the previous 2D case. 

Now, we have to rewrite and solve these equations in Eu- 
clidean space. We first note that there exists an analytic con- 
tinuation to Euclidean space yielding a finite instanton so- 
lution only for an open metric ansatz (k = —1) jsti- The 
equations of motion have the same form as equations (|6]l- 
(O after the limit T — > is taken, but with the substitution 
PCasA(3) -PcasA(3)- In Order to Calculate the instanton, 
we take the solution sitting at the unstable de Sitter vacuum, 
i.e., at the top of the /5casA(3) de Sitter hill (at R = Rq). The 
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equations now yield, 
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Since V(3)(i?o) > 0, these equations have a solution around 
the de Sitter vacuum given by. 
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withcj^ = pcasA(3)(-Ro)/(27ri?o^^pi)- The Hawking-Moss- 
like instanton can now be approximated by lUt], 



HM inst 
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Using i?o ~ 10" GeV"^ and V(3)(i?o) ^ GeV^ we 

obtain 5'HMinst — —10^^. The absolute value of this num- 
ber is much larger than the 2D case result from the previous 
section, however, it is still orders of magnitude smaller than 
—SdSi- The transition rate is again extremely suppressed. 

We will now search for finite interpolating solutions. Con- 
trary to the 2D case, most of the solutions are singular and 
one has to adopt a different strategy. We look for finite solu- 
tions around the minimum of the Euclidean potential, as was 
presented in Jst). From equation (|20] |. one obtains, 

R+2-R= (ipCa.A(3) - i^PcasA(3)) ■ (49) 



a 



After plugging R = R^ + SR and using relation ( |47] l. equation 
( |49] l becomes. 



SR + 2u!q cot(a;or) SR 
1 



PCasA(3) (^o) - |^0PcasA(3) (^o) 



SR = 0. (50) 



This equation has a solution expressed in terms of Gegenbauer 
polynomials iH, 



SR = A'Ci^) (cos(a;or)) , 



(51) 



where K <^ Rq is a constant, and the index satisfies the 
relation. 
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In order to have non-singular solutions, the Gegenbauer index 
n must be a non-negative integer. We find that this condition 
can only be satisfied for normal hierarchy Dirac neutrinos with 
the lightest neutrino mass ~ 7.4 x 10"^^ GcV. A quick check 
shows that this mass is associated with a stable 3D de Sitter 
vacuum, so the tunneling process is forbidden. 



However, as mentioned in fs'], it might be possible to obtain 
non-singular solutions in another way, i.e., by carefully fine- 
tuning the initial conditions rather than looking for a solution 
around the unstable de Sitter vacuum. Nevertheless, even if a 
non-singular solution exists, it is expected to yield an instan- 
ton action of the same order as the Hawking-Moss-like in- 
stanton action. Consequently, the transition rate is extremely 
small. 



IV. CONCLUSIONS 

In this paper, we have investigated the structure of the stan- 
dard model coupled to gravity at finite temperature for one- 
and two-dimensional compactifications (on a circle and on a 
torus). For each case, we have calculated the energy density 
and pressures from the free energy. Those quantities enter 
Einstein's equations, which govern the evolution of the com- 
pactified spacetime at finite temperature. We have found that 
stationary points of the equation of motion for the size of the 
compact space disappear with increasing temperature. The 
precise temperature at which this happens depends on the neu- 
trino masses, their hierarchy, and whether the neutrinos are 
Dirac or Majorana, but the qualitative behavior is the same. 

For the compactification on a circle we have shown numeri- 
cally, as an example, that the stable stationary point disappears 
at temperatures on the order of tens of Kelvin, for the lightest 
neutrino mass of ^ 10"^^ GeV, choosing the neutrinos to be 
Dirac with normal hierarchy. With increasing temperature, the 
stationary point never reappears. We have shown analytically 
that such behavior is expected in the high temperature limit. 
The case of the 2D compactification on a torus is very similar 

Finally, we have calculated the transition rates for tunnel- 
ing between vacua of different dimensionality at zero temper- 
ature. Following the steps outlined in jsl], we found Coleman- 
de Luccia-like instanton solutions for a tunneling from a 4D 
de Sitter universe to a spacetime with two spatial dimensions 
compactified on a torus. The rate for such a process turned 
out to be extremely suppressed. For the case of tunneling to 
a spacetime with one spatial dimension compactified on a cir- 
cle, we found only the Hawking-Moss-like instanton solution, 
which also yields a negligible transition rate. 

Nevertheless, under the condition that the universe has the 
right topology, because the transition rates at low tempera- 
tures are not zero, a tunneling process to a lower-dimensional 
spacetime might be the ultimate fate of our universe. 
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